Dark Matter annihilation in Draco: new considerations of the expected gamma flux 
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A new revision of the gamma flux that we expect to detect in Imaging Atmospheric Cherenkov 
Telescopes (IACTs) from neutralino annihilation in the Draco dSph is presented in the context of the 
minimal supersymmetric standard models (MSSM) compatible with the present phenomenological 
and cosmological constraints, and using the dark matter (DM) density profiles compatible with 
the latest observations. This revision takes also into account the important effect of the Point 
Spread Function (PSF) of the telescope, and is valid not only for Draco but also for any other DM 
target. We show that this effect is crucial in the way we will observe and interpret a possible signal 
detection. Finally, we discuss the prospects to detect a possible gamma signal from Draco for current 
or planned 7-ray experiments, i.e. MAGIC, GLAST and GAW. Even with the large astrophysical 
and particle physics uncertainties we find that the chances to detect a neutralino annihilation signal 
in Draco seem to be very scarce for current experiments. However, the prospects for future IACTs 
with upgraded performances (especially lower threshold energies and higher sensitivities) such as 
those offered by the CTA project, might be substantially better. 



I. INTRODUCTION 

Nowadays, it is generally believed that only a small 
fraction of the matter in the Universe is luminous. In the 
Cold Dark Matter (CDM) cosmological scenario around 
one third of the "dark" side of the Universe is supposed 
to be composed of weakly interacting massive particles 
(WIMPs). although other possible candidates like axi- 
nos or gravitinos are not excluded (see [1| for a recent 
review). The Standard Model can not provide a suitable 
explanation to the dark matter (DM) problem. However, 
its supersymmetric extension (SUSY) provides a natural 
candidate for DM in the form of a stable uncharged Ma- 
jorana fermion, called neutralino, which constitutes also 
one of the most suitable candidates according to the cur- 
rent cosmological constraints. At present, large effort is 
being carried out to detect this SUSY DM by different 
methods 0] ■ In the case of the new Imaging Atmospheric 
Cherenkov Telescopes (IACTs), the searches are based on 
the detectability of gamma rays coming from the anni- 
hilation of the SUSY DM particles. IACTs in operation 
like MAGIC or HESS Q, or satellites-based experi- 
ments like the upcoming GLAST satellite [5| , will play a 
very important role in these DM searches. 

A relevant question concerning the search of SUSY DM 
is where to search for the annihilation gamma ray signal. 
Due to the fact that the gamma flux is proportional to 
the square of the DM density, we will need to point the 
telescope to places where we expect to find a high con- 
centration of dark matter. In principle, the best option 
seems to be the Galactic Centre (GC), since it satisfies 
this condition and it is also very near compared to other 
potential targets. However, the GC is a very crowded 
region, which makes it difficult to discriminate between 
a possible 7-ray signal due to neutralino annihilation and 
other astrophysical sources. Whipple [f|, Cangaroo 0], 
and specially HESS [| and MAGIC [| have already car- 
ried out detailed observations of the GC and all of them 



reported a gamma point-like source at the Sag A* lo- 
cation. However, if this signal was interpreted as fully 
due to DM annihilation, it would correspond to a very 
massive neutralino very difficult to fit within the WMAP 
cosmology [l(| in the preferred SUSY framework [ll[ (al- 
though an alternative scenario with multi-TeV neutrali- 
nos compatible with WMAP is still possible, see [12j|). 
Furthermore, an extended emission was also discovered 
in the GC area, but it correlates very well with already 
known dense molecular clouds [l3|. Recently, new HESS 
data on the GC have been published and a reanalysis has 
been carried out by the HESS collaboration [l4j]. In this 
work, they especially explore the possibility that some 
portion of the detected signal is due to neutralino anni- 
hilation. According to their results, at the moment it is 
not possible to exclude a DM component hidden under 
a non-DM power-law spectrum due to an astrophysical 
source. 

There are also other possible targets with high dark 
matter density in relative proximity from us, e.g. the 
Andromeda galaxy, the dwarf spheroidal (dSph) galax- 
ies - most of them satellites of the Milky Way- or even 
massive clusters of galaxies (e.g. Virgo). DSph galaxies 
represent a good option, since they are not plagued by 
the problems of the GC, they are dark matter dominated 
systems with very high mass to light ratios, and at least 
six of them are nearer than 100 kpc from the GC (Draco, 
LMC, SMC, CMa, UMi and Sagittarius). 

Concerning DM detection, there are two unequivocal 
signatures that make sure that the 7-ray signal is due 
to neutralino annihilation: the spectrum of the source 
and its spatial extension. Certainly, the spectrum shape 
is crucial to state that the gamma source is due to DM 
annihilation, not only to discriminate it from the astro- 
physical backgrounds (e.g. hadronic, electronic and dif- 
fuse for IACTs), but also from any other astrophysical 
sources. The keypoint here is that the spectrum of a 
DM annihilation source is supposed to be in concordance 



2 



with that expected from the models of particle physics. 
Indeed, the exact annihilation spectra of MSSM neutrali- 
nos depend on the gaugino/higgsino mixing, but all of 
them show a continuum curved spectrum up to the mass 
of the DM particle and possibly faint 7-lines superim- 
posed [l5j . This spectrum will be very different from 
those measured in the case of an astrophysical source 
[li]]. As for the spatial extension of the source, in this 
case the important feature is that it should be extended 
and diffuse, showing also a characteristic shape of the 
gamma flux profile. Nevertheless, we must note that if 
we use an instrument that does not have a spatial res- 
olution good enough compared to the extension of the 
source, we might see only a point-like source instead of 
a diffuse or extended one. This means that, although 
we might reach high enough sensitivity for a successful 
detection, we would not be completely sure whether our 
signal is due to DM annihilation or not. Therefore, it is 
clear that it is really important to resolve the source so 
we can conclude that it can be interpreted as neutralino 
annihilation. 

In this work, we first calculate the gamma ray flux pro- 
files expected in a typical IACT due to neutralino anni- 
hilation in the Draco dSph, which represents in principle 
a good candidate to search for DM. Draco, located at 80 
kpc, is one of the dwarfs with many observational con- 
straints, which has helped to determine better its DM 
density profile. This fact is very important if we really 
want to make a realistic prediction of the expected 7-ray 
flux. These flux predictions have been already done for 
Draco u sing different models for the DM density profiles 
[3 EE Einiil- Nevertheless, in our case, we compute 
these flux predictions for a cusp and a core DM density 
profiles built from the latest stellar kinematic observa- 
tions together with a rigorous method of removal of in- 
terloper stars. This computation represents by itself a 
recommendable update of the best DM model for Draco, 
but as we will see it will be also useful to extract some im- 
portant conclusions concerning the possible uncertainties 
in the absolute 7-ray flux coming from astrophysics. 

Once we have obtained the flux profiles, we will use 
them to stress the role of the Point Spread Function 
(PSF) of the telescope. Including the PSF, which is di- 
rectly related to the angular resolution of the IACT, is 
essential to interpret correctly a possible signal profile 
due to neutralino annihilation, not only for Draco but 
also for any other target. In fact, we will show that, de- 
pending on the PSF of the IACT, we could distinguish or 
not between different models of the DM density profile 
using the observed flux profile. In the case of the cusp 
and core DM density profiles that we use, it could be im- 
possible to discriminate between them if the PSF is not 
good enough. It is worth mentionin g th at most of previ- 
ous works in the literature (except [2l|) that calculated 
the expected flux profiles in IACTs due to dark matter 
annihilation did not take into account this important ef- 
fect. Because of that, to emphasize the role of the PSF 
constitutes also one of the goals of this work. 



Finally, we present the DM detection prospects of 
Draco for some current or planned experiments, i.e. 
MAGIC, GLAST and GAW. We carry out the calcula- 
tions under two different approaches: detection of the 
gamma ray flux profile from the cusp and core DM mod- 
els for Draco, and detection of an excess signal in the 
direction of the dwarf galaxy. We will show how the first 
approach could give us a lot of information about the 
origin of the gamma ray flux profile, but it is slightly 
harder to have success for this case compared to the sec- 
ond approach, where even the PSF of the instrument is 
not essential and still we could extract some important 
conclusions if we reached the required sensitivity. 

The paper is organized as follows. In Section|TI]we first 
present all the equations necessary to properly calculate 
the expected 7-ray flux in IACTs due to neutralino anni- 
hilation. Both the particle physics and the astrophysics 
involved are described carefully. For the particle physics 
we analyze neutralino annihilation in the context of the 
minimal supersymmetric standard models (MSSM) com- 
patible with the present phenomenological and cosmo- 
logical constraints. In Section ITU we show in detail the 
model that we use for the DM distribution in Draco. In 
Section ITVl we calculate the Draco flux predictions. We 
also stress the important role of the PSF. In Section fVl 
the prospects to detect signal due to neutralino annihi- 
lation in Draco are shown for some current or planned 
experiments, i.e. MAGIC, GLAST and GAW. Conclu- 
sions are finally given in Section IVIl 

II. THE 7-RAY FLUX IN IACTs 

The expected total number of continuum 7-ray pho- 
tons received per unit time and per unit area in a tele- 
scope with an energy threshold E t h is given by the prod- 
uct of two factors: 

F(E > E th ) = ^-fsusr ■ U(* ). (1) 

where represents the direction of observation rela- 
tive to the centre of the dark matter halo. The factor 
fsuSY includes all the particle physics, whereas the fac- 
tor U(9 ) involves all the astrophysical properties (such 
as the dark matter distribution and geometry consider- 
ations) and also accounts for the beam smearing of the 
telescope. 

A. Particle physics: the fsuSY parameter 

In R-parity conserving supersymmetric theories, the 
lightest SUSY particle (LSP) remains stable. The widely 
studied Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM) can predict a neutralino as the LSP 
with a relic density compatible with the WMAP bounds. 
In this section we concentrate on the MSSM under the 
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assumption that the neutralino is the main component 
of the DM present in the universe. We consider its abun- 
dance inside the bounds 0.09 < Q x h 2 < 0.13 as derived 
by fitting the ACDM model to the WMAP data [ijj, 
although the possibility for other DM candidates is not 
excluded (see Ref . [3, S| and references therein) . 

The properties of the neutralinos in the MSSM are 
determined by its gaugino-higgsino composition: 

X = Xi = N n B + N 12 W 3 + N 13 H° + N 13 H% (2) 

At leading order, neutralinos do not annihilate into two- 
body final states containing photons. However at one 
loop it is possible to get processes such as 0, SS S3 : 

X + X -> 77 
X + X -> Zj 

with monochromatic outgoing photos of energies 

o 

Tfl 

E~ ~ m v , E~ <~~i m v — (3) 

4m x 

respectively. Also, the neutralino annihilation can pro- 
duce a continuum 7 - ray spectrum from hadronization 
and subsequent pion decay which can be dominate over 
the monochromatic 7's. 

The lack of experimental evidence of supersymmetric 
particles leaves us with a number of undetermined pa- 
rameters in the SUSY models. Therefore, the annihila- 
tion cross sections involved in both the computation of 
Q x h 2 and 7 production can change orders of magnitude 
with neutralino mass and its gaugino-higgsino composi- 
tion. To be specific we consider mSUGRA models, where 
the soft terms of the MSSM are taken to be universal at 
the gauge unification scale Mqut- Under this assump- 
tion, the effective theory at energies below Mqut de- 
pends on four parameters: the soft scalar mass mo, the 
soft gaugino mass mi/2, the soft trilinear coupling Aq, 
and the ratio of the Higgs vacuum expectation values, 
tan/3 = (H®) / (H%). In addition, the minimisation of 
the Higgs potential leaves undetermined the sign of the 
Higgs mass parameter /i. 

To provide some specific values, we assume Aq = 0, 
/i > and two values of tan/3, 10 and 50. With this 
two cases we can provide a qualitative picture of most 
relevant space of parameters and the constraints imposed 
by phenomenology. The impact of the size of Ag or the 
sign of /i can be found in several studies of the parameter 
space of the MSSM (see for example Refs. [M iM S3, HI 
and references therein). 

In Fig. Q] we displayed some lines of constant values of 
2 (i><7 77 ), (va-yz) and cross sections for continuum 7— ray 
emission on the plane mo - along with the con- 

straints derived from the lower bound on the mass of the 
lightest neutral Higgs, m^ — 114.1 GeV, chargino mass 
m x + = 103 GeV and BR(6 — > 37), and the areas with 
Qh 2 on the WMAP bounds. Also we provide the lines 
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FIG. 1: Contours on the mo — rax/i plane, the up and down 
ruled areas are excluded by the not satisfaction of the EWSB 
(up) and because m x > mf. The area below the upper thick 
solid line satisfies the experimental bound on the chargino 
mass, while the green shaded areas indicate the areas that 
predict neutralino relic density on WMAP bounds. From left 
to right, the ruled areas are excluded by the bounds on m° and 
BR(b — > S7) respectively. The doted lines indicates the values 
of a xp in pb, the dash and dot-dash lines corresponds respec- 
tively to 2va xx —,~ n and va xx ^z-y in units of 10 -29 cm 3 s _1 . 
The line of crosses (plus) correspond to the continuum pho- 
ton production with E-y > 1 GeV (E-y > 100 GeV) also in 
units of 10~ 29 cm 3 s~\ 

with constant values for the elastic scattering %-proton, 
relevant for neutralino direct detection. In the computa- 
tion we used DarkSUSY 29] combined with isasugra [3(| 
implementing the phenomenological constraints as dis- 
cussed in Ref. [26[; the estimation of BR(6 — » sj) was 
performed using Refs. [3lT I32I] . 

On the lower area consistent with WMAP the neu- 
tralino is Bino-like, the relic density is satisfied mostly 
due to coannihilations x — t and in the case of tan /3 = 50 
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because of annihilation \ — \ through resonant chan- 
nels. In this sector (2vct 77 } is dominant by a factor of 10 
with respect to (va~ t z), however its larger values lie in 



the areas constrained by the bounds on to 



X+i m h° 



and 



BR(b -> s 7 ). 

The higher area consistent with WMAP lies on the 
hyperbolic branch, the neutralino is gaugino-higgsino 
mixed. The position of this region is very dependent 
on the mass of the top; we used m t = 175 GeV. 

In Fig. [5] we present the values of fsusY versus the 
threshold energy of the detector, starting at E th — 
0.4 GeV, for neutralinos satisfying relic density and phc- 
nomenological constraints. fsuSY is calculated as: 



SUSY 



0(E th > m x ) ■ 2 (wcr 77 ) 



2m\ 



9{E t h > m x - 3^f-) • (va lZ ) + k (va cont .) 



2m x 



(4) 



where 9 is the step function and k the photon multiplicity 
for each neutralino annihilation. We display in different 
panels values of fsuSY corresponding to points on hyper- 
bolic branch (mo > 2 TeV as we can see in FigJT]) and the 
values corresponding to the x — t coannihilation and reso- 
nant annihilation areas (mo < 2 TeV) . We can appreciate 
that on the x ~ t coannihilation area, m x has an upper 
bound beyond which the relic density constraint is no 
longer satisfied while on the hyperbolic branch no upper 
bound for m x is reached. The largest values of fsuSY are 
not present since they correspond to low values of m 1 / 2 
suppressed by the bounds on m x + , m h o and b — ► sj. It 
is interesting to remark that the higher values of fsuSY 
on the constrained areas lie in which the <r x _ p reaches 
values in the range of direct detection experiments like 
Genius 

To show a more general sampling on the parameter 
space, we present in Fig. [3] the 7-ray production in neu- 
tralino annihilation versus m x . Each panel contains 
about 1400 points satisfying the WMAP and phenomeno- 
logical constraints mentioned above. These points are se- 
lected from random scan over 600000 parameter on the 
ranges: 



3 < tan/3 < 60, 

50 GeV < mi/2 < 2300 GeV, 

—3 too < Aq < 3 mo, 

50 GeV < too < 10 TeV, (5) 

and both signs of the [i term. 

In the top panel of Fig. [3] we include models with 
m < 2 TeV. These points satisfy the relic density con- 
straint due mostly to f — x coannihilations, and to a 
minor extend due to resonances in the annihilation chan- 
nels (for large values of tan /3); m x remain below 700 GeV 
in this panel. The continuum emission dominates for 



Ej > 1 GeV, while for E y > 100 GeV we find points 
where the continuum production is of the same order as 
the monochromatic 7's. The bottom panel of Fig. [3]con- 
tains models too > 2 TeV. These correspond mostly to 
the hyperbolic branch region; also we can appreciate that 
larger values for m x are allowed in this area than in the 
top panel. The points with larger cross sections in both 
panels corresponds to the smaller values of the pseudo- 
scalar Higgs mass ttla and higher higgsino composition of 
the neutralino. These two factors favour the fermion pro- 
duction in neutralino annihilation channels. On the top 
panel, the larger cross sections correspond to the larger 
values of tan j3 while on the right panel low values for to a 
can be reached for any tan [3. 

The mSUGRA results do not differ significantly that 
the ones obtained in the more general MSSM models ob- 
tained by waiving the universality conditions on the soft 
terms as we will present in section I"VB1 



B. Astrophysics: the U^o) parameter 

All the astrophysical considerations are included in the 
expression U(if>o) in Eq.JT]). This factor accounts for the 
dark matter distribution, the geometry of the problem 
and also the beam smearing of the IACT, i.e. 



U(V ) = J J(^)B(Q)dQ 



(6) 



where B(il)dfl represents the beam smearing of the tele- 
scope, commonly known as the Point Spread Function 
(PSF). The PSF can be well approximated by a Gaus- 
sian: 



B(n)dfl = exp 



2o1 



sin d9 



with at the angular resolution of the IACT. It is worth 
mentioning that there is some dependence of a t with the 
energy (see e.g. [13]). However, for simplicity, we will 
suppose this term to be constant. The PSF plays a very 
important role in the way we will observe a possible DM 
signal in the telescope. However, most of previous works 
in the literature did not take into account its effect (ex- 
cept [2l|]; in [11] the PSF apparently was also used, al- 
though it is not mentioned in the text (S. Profumo, pri- 
vate communication)). In Section lIVI we will study in 
detail the importance of the PSF in the determination of 
the gamma ray flux profile. 

The J (if?) factor of Eq.© represents the integral of 
the line-of-sight of the square of the dark matter density 
along the direction of observation if?: 



J(if?) = 



2 

Pdr, 



(r) dX 



2 

Pdr, 



[r(X)} dX (8) 
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FIG. 2: Values of fsuSY respect Eth, for the points in Fig[T]on the WMAP region and satisfying all the phenomenological 
constraints. The ruled areas include the continuum 7— ray emission with E-, > Eth, while the ones limited by thin lines 
correspond only to monochromatic channels. The left (right) panel corresponds to the upper (lower) allowed WMAP area in 
Fig. 1. 



Here, r represents the galactocentric distance, related to 
the distance A to the Earth by: 



r = J\ 2 + R% - 2 A R Q cos * (9) 



where Rq is the distance from the Earth to the centre of 
the galactic halo, and \1/ is related to the angles 9 and (j> by 
the relation cos ^ = cos ipo cos 9 + sin t/jg sin 9 cos <fi. The 
lower and upper limits X m in and \ ma x in the l.o.s. inte- 
gration are given by Rq cos ij) ± \J r 2 — Rq sin 2 ip, where 
rt is the tidal radius of the dSph galaxy in this case. 



III. DARK MATTER DISTRIBUTION IN 
DRACO 

In our modelling of Draco we used the sample of 207 
Draco stars with measured line-of-sight velocities origi- 
nally considered as members by (34[- In selecting these 
stars these authors relied on a simple prescription going 
back to [HI and based on rejection of stars with veloci- 
ties exceeding 3<ti os where o\ QS is the line-of-sight velocity 
dispersion of the sample. [36| have shown that if all these 
207 stars are used to model Draco velocity distribution 
the resulting velocity moments can be reproduced only 
by extremely extended mass distribution with total mass 
of the order of a normal galaxy. Their arguments strongly 
suggested that some of the stars may in fact be unbound 
and the simple 3cri os rejection of stars is insufficient. 

Here we apply a rigorous method of removal of such 
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FIG. 3: Values of n 7 < a xx v > in err? j s including continuum 
emission for E-, > 1 GeV (circles), _E 7 > 100 GeV (stars), and 
considering only the two monochromatic channels (plus). 



interlopers originally proposed by [37J and applied to 
galaxy clusters. The method relies on calculating the 
maximum velocity available to the members of the ob- 
ject assuming that they are on circular orbits or infalling 
into the structure. The method was shown to be the most 
efficient among many methods of interloper removal re- 
cently tested on cluster-size simulated dark matter haloes 
by [38|] . Its applicability and efficiency in the case of dSph 
galaxies was demonstrated by (39|. Fig. [4] shows the re- 
sults of the application of this procedure to Draco. The 
207 stars shown in the plot are divided into those itera- 
tively rejected by the procedure (open circles) and those 
accepted at the final iteration (filled circles). 

The final sample with 194 stars is different from any 
of the three considered by [36| therefore we repeat their 
analysis here for this new selection. Our analysis is ex- 
actly the same, except that in the calculation of the veloc- 
ity moments we use 32-33 stars per bin instead of about 
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FIG. 4: The line-of-sight velocities versus projected distances 
from the galaxy centre for 207 stars from [3J|. Open circles 
mark the 13 stars rejected by our interloper removal proce- 
dure, filled circles show the 194 ones accepted. 



40 and we consider a DM profile with a core in addition 
to the cusp one. The profiles of the line-of-sight velocity 
moments, dispersion and kurtosis, obtained for the new 
sample are shown in Fig. [5l The kurtosis was expressed 
in terms of the variable k = [log(3if/2.7)] 1 / 10 where K 
is the standard kurtosis estimator. We assumed that the 
DM distribution in Draco can be approximated by 



p d (r) = Cr "exp 



r 



(10) 



proposed by [40j, which was found to fit the density dis- 
tribution of a simulated dwarf dark matter halo stripped 
during its evolution in the potential of a giant galaxy. In 
the same work it was found that the halo, which initially 
had a NFW distribution, preserves the cusp in the inner 
part (so that a = 1 fits the final remnant very well) but 
develops an exponential cut-off in the outer parts. Here 
we will consider two cases, the profile with a cusp a = 1 
and a core a = 0. It remains to be investigated which 
scenarios could lead to such core profiles. 

The best-fitting solutions to the Jeans equations (see 
[36| ) for two component models with dark matter profiles 
given by (fTU|) are plotted in Fig. \5\ as solid lines in the 
case of the cusp profile and dashed lines for the core. 
The best-fitting parameters of the two models are listed 
in Table U where Md /Ms is the ratio of the total dark 
matter mass to total stellar mass, r^/Rs is the break 
radius of equation (fTTJj) in units of the Sersic radius of 
the stars and is the anisotropy parameter of the stellar 
orbits. 

Fig. [S] shows the best-fitting dark matter density pro- 
files in the case of the cusp (solid line) and the core 
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FIG. 6: The best-fitting DM density profiles for Draco with 
a cusp (solid line) and a core (dashed line). 



(dashed line). As we can see, both density profiles are 
similar up to about 1 kpc, where they are constrained 
by the data. The reason for very different values of the 
break radius in both cases is the following. The kur- 
tosis is sensitive mainly to anisotropy and it forces f3 to 
be close to zero in both cases. However, to reproduce the 
velocity dispersion profile with /3 ~ the density profile 
has to be steep enough. In the case of the core it means 
that the exponential cut-off has to occur for rather low 
radii, which is what we see in the fit. The cusp profile 
does not need to steepen the profile so much so it is much 
more extended and its total mass is much larger. 



IV. DRACO GAMMA RAY FLUX PROFILES 



FIG. 5: The line-of-sight velocity dispersion (upper panel) 
and kurtosis variable k (lower panel) calculated for the sample 
of 194 stars with 32-33 stars per bin. The lines show the best- 
fitting solutions of the Jeans equations for the DM profile with 
a cusp (solid lines) and a core (dashed lines). 



TABLE I: Best-fitting parameters of the two-component mod- 
els for the DM profiles with a cusp (a = 1) and a core (a = 0) 
obtained from joint fitting of velocity dispersion and kurtosis 
profiles shown in Fig. \S\ The last column gives the goodness 
of fit measure x 1 /N. 



profile 


M D /Ms 


r h /Rs 





X 2 /N 


cusp 


830 


7.0 


-0.1 


8.8/9 


core 


185 


1.4 


0.06 


9.5/9 



In order to compute the expected gamma flux, we 
need to calculate the value of the "astrophysical factor" , 
U(ipo), given in Eq.{T]) and presented in detail in Section 
III Bl We calculated it for the core (a — 0) and cusp 
(a = 1) density profiles given by Eq. fTU)) using the pa- 
rameters listed in Table HT1 (that were deduced from those 
given in Tableland where we used Rs = 7.3 arcmin for 
Draco, following [ll|). i? Q was set to 80 kpc, as derived 
from an analysis on the basis of wide-field CCD photom- 
etry of resolved stars in Draco [52]. For the tidal radius 
we used a value of 7 kpc as given by [l6| and derived 
from the Roche criterion supposing an isothermal pro- 
file for the Milky Way. Nevertheless, this value depends 
strongly on the profile used for the Milky Way and Draco, 
e.g. a value of 1.6 kpc is found when a NFW DM density 
profile is used for both galaxies [lfij . It is worth men- 
tioning, however, that we computed J(^E') for different 
values of r t and we checked that the difference between 
choosing r t = 1.6 kpc and rt = 7 kpc is less than 5% for 



TABLE II: Values of C and rj, for a cusp and a core DM 
density profile given by Eq.fTUJ), as deduced from those pa- 
rameters listed in Table |U 



profile 




C 


r b (kpc) 


cusp 


3.1 x 


10 Y Mo/kpc^ 


1.189 


core 


3.6 x 


10 8 M Q /kpc 3 


0.238 



* = 0.5°, and still less than 10% for * = 1°. Therefore, 
for the case of Draco, any value r t > 1.5 leads to robust 
and very similar results. 

There is another issue that we will have to take into 
account in order to calculate U(ipo)- If we integrate the 
square DM density along the line of sight using a cusp 
DM density profile, we will obtain divergences at angles 
i/jq — » (clearly there will not be any problem for core 
profiles). This can be solved by introducing a small con- 
stant DM core in the very centre of the DM halo. In 
particular, the radius r cut at which the self annihilation 
rate ti ~ (< a ann v > n x r cu t) _1 equals the dynamical 
time of the halo td yn ~ (G p) -1 / 2 , where p is the mean 
halo density and n x is the neutralino number density, is 
usually taken as the radius of this constant density core 
[3 . For the NFW DM density profile this value for r cut 
is of the order of 10~ 13 — 10~ 14 kpc. For steeper DM den- 
sity profiles (such as the compressed NFW or the Moore 
profile) a value of r cut ~ 10 -8 kpc is obtained. We used 
a value of 10~ 8 kpc in all our computations. We must 
note that r cu t represents a lower limit concerning the ac- 
ceptable values for this parameter, so the obtained fluxes 
should be taken as upper bounds. 

Once we have calculated U(ipo), we will need also 
to take a value for the fsusY parameter in order 
to obtain the absolute flux due to neutralino annihi- 
lation (see Eq. [!}. We chose a value of fsusY = 
1CP 33 ph GeV~ 2 cm 3 s _1 in all our computations for 
a typical E th ~ 100 GeV of the IACT. In the framework 
of MSSM, this value corresponds to one of the most op- 
timistic values that is possible to adopt for fsusY ac- 
cording to Fig. [2] for the two different values of tan (3 
presented. 

The resulting 7-ray flux profiles for Draco are plotted 
in Figure [3 where we used a PSF with a t = 0.1° (to 
simplify the notation, hereafter we will use PSF= 0.1° 
to refer to a PSF with a t = 0.1°). This value of 0.1° 
is the typical value for an IACT like MAGIC or HESS. 
It is important to note that the core and a cusp density 
profiles would be distinguishable thanks to a different 
and characteristic shape of the flux profile in each case. 

To illustrate the PSF effect on the shape of the ob- 
served flux profile with IACTs, in the top panel of Figure 
|5]we show the same as in Fig. but here for a PSF= 1°. 
As we can see, although we have different DM density 
profiles, a worse telescope resolution makes both result- 
ing flux profiles for a core and a cusp indistinguishable. 
We may think that we could distinguish them from the 
value of the absolute flux. However, the difference in the 
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FIG. 7: Draco flux predictions for the core (dashed line) and 
cusp (solid line) density profiles, computed for a typical IACT 
with E th = 100 GeV and a PSF= 0.1°. We used fsusv = 
10~ 33 ph GeV~ 2 cm 3 s" 1 (see text for details). 



absolute flux between both DM density profiles is very 
small and in practice the distinction would be impossi- 
ble. There are many uncertainties in this absolute flux 
coming from the particle physics. fsusY may be very dif- 
ferent from the most optimistic case assumed here, since 
it could vary more than three orders of magnitude for 
this SUSY model (see Fig. |5J). On the other hand, the 
uncertainty in flux due to the DM density profile to be 
core or cusp is negligible at least in the inner 0.5 degrees. 

Concerning the effect of the PSF given the same DM 
density profile, a worse telescope resolution flattens the 
flux profile as expected. It can be clearly seen in the 
bottom panel of Figure [51 where we plot the Draco 7-ray 
flux predictions only for the cusp density profile but using 
two different values of the PSF (0.1° and 1°), and where 
we plot also the same flux profile computed without PSF 
for comparison. 

A good example to show the real importance of the 
telescope resolution can be found around the controversy 
generated in the wake of the Draco 7-ra y ex cess reported 
by the CACTUS collaboration in 2005 At this mo- 
ment, it seems clear that this excess was not real and was 
probably due to instrumental and trigger- related issues. 
However, concerning our line of work and always just 
with the intention of clarifying the role of the PSF, we 
must mention the results shown in [3- There, the CAC- 
TUS data were superimposed on different flux profiles 
(each of them related to possible models of DM density 
profiles for Draco) in Figure 2. These flux profiles were 
computed using an angular resolution of 0.1°, whereas 
the CACTUS PSF is quite worse than that (around 0.3° 
for the Crab and probably worse for Draco [H, [3)- 
Looking at that figure (despite the authors' indication 
to take it with care) one may come to the conclusion 
that a core profile seems to be the most adequate DM 
density profile for Draco, as opposed to the cusp profile. 
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FIG. 8: Top panel: Draco flux predictions for the core (dashed 
line) and cusp (solid line) density profiles, computed using 
a PSF= 1°. Bottom panel: Draco flux predictions for the 
cusp density profile using two different PSFs. Solid line cor- 
responds to PSF=0.1° and dashed line to PSF=1°. The flux 
profile computing without PSF is also shown for comparison 
(dot-dashed line). 



However, it would be more appropriate to make this com- 
parison between the CACTUS data and the flux profiles 
using in both cases the same PSF of the experiment. Do- 
ing so and taking into account the PSF effect properly, it 
would be difficult to use the CACTUS data to discrimi- 
nate between different models for the DM density profile 
as described in [33| , since all of the resultant flux profiles 
would have essentially the same shape. Only the abso- 
lute flux could give us a clue to make this distinction 
possible, but as already mentioned there are too many 
uncertainties related to an absolute value to be able to 
extract solid conclusions. 



V. DETECTION PROSPECTS FOR SOME 
CURRENT OR PLANNED EXPERIMENTS 

A. Flux profile detection 

Would it be possible to detect a signal due to neu- 
tralino annihilation in Draco using present or planned 
IACTs and satellite-based gamma ray experiments? Al- 
though there are many uncertainties concerning both 
particle physics and astrophysical issues, as already 
pointed out, it is possible (and necessary) to make some 
calculations. These calculations will allow us to estimate 
at least the order of magnitude of the flux that we could 
expect in our telescopes, and will help us learn which 
instrument is, in principle, the best positioned and opti- 
mised to detect a possible signal from Draco. 

Draco is located in the northern hemisphere, more pre- 
cisely at declination 58°. Because of that, and regard- 
ing currently operating IACTs, MAGIC 3 and VERI- 
TAS ____[ are the best options thanks to their geograph- 
ical position. Since both experiments have comparable 
sensitivity and PSF, we will focus only on MAGIC. For 
this IACT, with an energy threshold around 50 GeV for 
zenith observations, Draco could be observed 61° above 
the horizon, so an energy threshold ~ 100 GeV seems to 
be still possible at that altitude. 

Concerning GLAST [f| , this satellite-based experiment 
is designed for making observations of celestial gamma- 
ray sources in the energy band extending from ~ 10 
MeV to 300 GeV, which is complementary to the one for 
MAGIC. Moreover, it will have a PSF ~ 0.1° at 10 GeV 
[47| . which will make this instrument very competitive 
also in DM searches. GLAST is planned to be launched 
at the end of 2007. 

The flux profiles detection prospects for both gamma- 
ray experiments can be seen in Figure ___ where the lines 
of sensitivity for MAGIC and GLAST are superimposed 
on the flux profile computed using the cusp DM density 
profile for Draco as given by Eq. (fT0|) together with the 
parameters listed in Table HT1 For the case of MAGIC, 
the sensitivity line represents 250 hours of observation 
time and a 5<7 detection level. As pointed out before, 
this curve should be comparable and valid for VERITAS 
experiment as well. Concerning the predicted sensitivity 
for GLAST, it was calculated by the GLAST team 
for 1 year of observation and a 5cr detection level. Val- 
ues of fsusY = 10~ 33 ph GeV 2 cm 3 s" 1 at 100 GeV 
(MAGIC) and Jsusy = 1 • 10" 29 ph GeV" 2 cm 3 s" 1 at 1 
GeV (GLAST) were chosen to convert the values calcu- 
lated from the astrophysical factor U(^fo) to flux. These 
values of fsuSY correspond to the most optimistic values 
compatible with that shown in Fig. ___ We chose E t h = 1 
GeV for GLAST to be sure that we will have a good 
angular resolution for the arriving photons and also to 
avoid some uncertainties at low energies in the response 
of the GLAST LAT detector as pointed out in [11]. 

From Figure ___ we can see that the detection of the 
gamma ray flux profile due to neutralino annihilation 



10 




0.0 0.3 0.6 0.9 1.2 

fo (deg.) 




0.0 0.3 0.6 0.9 1.2 

fo (deg.) 



FIG. 9: Draco flux profile detection prospects for MAGIC 
(top panel) and GLAST (bottom panel). The flux profile 
(solid line) corresponds to the cusp density profile given in 
Table [TT] and using a PSF=0.1°. The sensitivity lines for both 
instruments were computed for a 5a detection level and 250 
hours (MAGIC) / 1 year (GLAST) of integration time. These 
figures should be taken as the most optimistic case, since we 
used the most optimistic scenario from particle physics in each 
case. We adopt E th = 100 GeV for MAGIC and E th = 1 GeV 
for GLAST. See text for details. 



in Draco with the MAGIC telescope seems to be im- 
possible (at least for the DM density profiles and the 
particle physics model used here), since we would need 
roughly five orders of magnitude more sensitivity than 
that reached by this instrument in a reasonable time to 
have at least one opportunity to detect a possible signal 
coming from DM annihilation. For the case of GLAST 
we obtain more or less the same, although in this case the 
expectations are substantially better and we would need 
three orders of magnitude more sensitivity. We must note 
that this is the most optimistic scenario, so even leaving 
space to take into account the large uncertainties coming 
from the particle physics it seems very hard to have any 
chance of detection. In fact, in case of adopting a pes- 
simistic value for fsuSY, the 7-ray flux profile shown in 



Fig. [9] could decrease in more than three orders of mag- 
nitude easily. 

There are some issues concerning a possible detection 
of DM annihilation not only in Draco but also in any 
other possible target that should be taken into account 
at this moment. In the case of a positive detection, this 
signal would be diffuse (i.e. no point-like source) for an 
instrument with a PSF good enough. This fact and the 
characteristic spectrum of the source represent the best 
clues to distinguish between a 7-ray signal due to neu- 
tralino annihilation from other astrophysical sources. In 
the case of Draco, however, we are very far from obtain- 
ing any of these clues, so Draco does not seem a good 
target for DM searches for current or planned 7-ray ex- 
periments (MAGIC-II, for example, will have a factor of 
2 more sensitivity than the single MAGIC telescope , 
still clearly insufficient, and the same improvement in 
sensitivity is reached in the case of the GLAST satellite 
if we observe 5 years instead of 1). 



B. Excess signal detection 

Although it is strongly recommendable to try to detect 
the gamma ray flux profile due to neutralino annihilation, 
we could also search only for a gamma ray excess signal 
in the direction of Draco with a DM origin. In this case 
the prospects should be somewhat better, since we are 
only interested in flux detectability (i.e. no interested 
in observing an extended emission or not; no discrimi- 
nation between different flux profiles). Draco is around 
2° in the sky (the whole DM halo), and one MAGIC 
pointing is 4°x 4°, so the whole galaxy is within one of 
these MAGIC pointings. This means that, to calculate 
the prospects of an excess signal detection for MAGIC, 
we must sum the gamma ray fluxes coming from all the 
neutralino annihilations that occur in the entire halo of 
the dSph. The same is valid for the case of GLAST, since 
one of the main objectives of this mission will be to sur- 
vey the whole sky, so Draco as a whole will be observed. 
Nevertheless, most of the 7-ray flux due to DM annihi- 
lation comes from the inner regions of the dwarf, so it 
would be better to integrate the flux and make all the 
calculations only for those regions. Otherwise, if we inte- 
grate up to very large angles from the center, we would 
be increasing the noise in a large amount without prac- 
tically increasing the 7-ray signal, since the gamma ray 
flux profile decreases very rapidly from the center to the 
outskirts. 

Wc would like here to take the opportunity to mention 
GAW [50j, which is a R&D path- finder experiment, still 
under development, for wide field 7-ray astronomy. GAW 
will operate above 0.7 TeV and will have a PSF ~ 0.2°. It 
will consist of three identical telescopes working in stereo- 
scopic mode (80m side). The main goal of GAW is to 
test the feasibility of a new generation of IACTs, which 
join high sensitivity with large field of view (24° x 24°). 
GAW is planned to be located at Calar Alto Observatory 
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TABLE III: Prospects of an excess signal detection for 
MAGIC, GLAST and GAW. Concerning the integrated flux 
for Draco, Foraco, the most optimistic and pessimistic values 
are given in the form Fn ra co,min - F Dra co, m ax- F min repre- 
sents the minimum detectable flux for each instrument. Both 
the integrated absolute fluxes for Draco and the values given 
for the sensitivities (i.e. the minimum detectable flux, F m i n ) 
refer to the inner 0.5° of the galaxy (although its size is ~ 
2°), just to improve the signal to noise ratio as explained in 
the text. 
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(Spain) and a first part of the array should be completed 
and start to operate during 2008. It is a collaborative 
effort of research Institutes in Italy, Portugal, and Spain. 
We will also present some calculations concerning the 
possibility to observe a 7-ray signal in the direction of 
Draco by GAW, just to illustrate the capabilities of the 
instrument. Nevertheless, the main advantages of GAW 
will point to other directions, e.g. the possibility to sur- 
vey a large portion of the sky in a reasonable time above 
0.7 TeV. 

In Table IIIII we show the prospects of an excess signal 
detection (5er level) for MAGIC, GLAST and GAW and 
for the case of the cusp DM density profile. Both the 
integrated absolute fluxes for Draco and the values given 
for the sensitivities (i.e. the minimum detectable flux, 
Fmin) refer to the inner 0.5° of the galaxy (although the 
total size of Draco in the sky is ~ 2°), just to improve 
the signal to noise ratio as explained above. The inte- 
grated absolute flux for Draco is not the same for the 
three experiments, since the fsusY parameter depends 
on the energy threshold of each instrument, which is dif- 
ferent (we chose 100 GeV for MAGIC, 1 GeV for GLAST 
and 700 GeV for GAW). We did the calculations for the 
most optimistic case (i.e. the highest value of fsuSY that 
we could adopt in the MSSM scenario following Fig. [21 
and given the energy threshold of each telescope) and the 
most pessimistic one (the lowest Jsusy)- 

According to Table IIIII an excess signal detection 
would be impossible with any of these three instruments, 
since even in the most optimistic cases they do not reach 
the required sensitivity by far. If we compare these num- 
bers with the maximum fluxes shown in Figj9l we can 
see that the enhancement, although relevant (roughly a 
factor of 4 in the best case) , is still very insufficient for a 
successful detection. 

Moreover, from the results shown in Table IIIII it be- 
comes really hard to extract any useful results that would 
help us understand at least a bit better the problem of the 
dark matter. Indeed, if we observe Draco with MAGIC 
or GLAST and we do not detect any excess signal, we 
will not be able to put any useful constraints on the par- 
ticle physics involved. None of the actual allowed values 



for fsuSY could be rejected if we find no detection. This 
makes Draco even less attractive for those current gamma 
ray experiments that try to find DM traces. 
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FIG. 10: Exclusion limits for MAGIC (upper panel) and 
GLAST (bottom panel), for continuum 7-ray emission above 
100 GeV (MAGIC) and 1 GeV (GLAST). The lines repre- 
sent for each instrument the minimum detectable n 7 < av > 
adopting the cusp DM density profile given in Table [TT] for 
Draco. We used 250 hours (1 year) of integration time for 
MAGIC (GLAST) and a 5<r detection level. Below the lines, 
the SUSY models (points) do not yield a flux high enough 
for a successful detection. According to these plots, neither 
MAGIC nor GLAST may put any useful constraints at least 
to this particular particle physics model. 

This can be clearly seen in Fig. [TQl where we show 
the parts of the SUSY parameter space that will be de- 
tectable by MAGIC and GLAST for the case of Draco 
and adopting the cusp DM density profile presented in 
Section IIIII The points represent MSSM models, in con- 
trast with Fig. [3] the mSUGRA condition of Eq. O has 
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been waived such that models with random sfermion 
masses below 10 TeV and gaugino masses below 3 TeV 
are included (we considered equal soft terms for the first 
two generations to avoid contradiction with flavor vio- 
lating observables). In the case of MAGIC we only plot 
those values of n 7 < av > computed for E 1 > 100 GeV 
(which is the MAGIC E th ), and for GLAST only those 
ones computed for E 1 > 1 GeV (GLAST E th ). We must 
note that for each case we include in the same figure all 
the points no matter the value of Too ( a distinction was 
done in Fig. [3]). Also, the points related to the monochro- 
matic channels are not shown in both figures, since their 
values are negligible comparing to those coming from the 
continuum emission and therefore they are not relevant 
here. 

The limit between the detectable and the non- 
detectable areas is given by the dashed lines, so that those 
SUSY points that lie above these lines yield a detectable 
flux and the points below them are not accessible to ob- 
servation. As already pointed above, it is clear that the 
constraints imposed by both experiments to this particu- 
lar particle physics model (MSSM) are very relaxed and 
in fact the detection lines do not reach in any case any of 
the SUSY models by more than two orders of magnitude 
in the best case. Again, as in Fig. [51 the prospects for 
GLAST are somewhat better than for MAGIC, but they 
are still clearly insufficient to extract relevant results or 
conclusions. 



VI. CONCLUSIONS 

In this work we focused on the possibility to detect a 
signal coming from neutralino annihilation in the Draco 
dwarf. This galaxy, a satellite of the Milky Way, rep- 
resents one of the best suitable candidates to search for 
dark matter outside our galaxy, since it is near and it has 
probably more observational constraints than any other 
known dark matter dominated system. This fact be- 
comes crucial when we want to make realistic predictions 
of the expected observed 7-ray flux due to neutralino an- 
nihilation. 

Draco is a dwarf galaxy tidally stripped by the Milky 
Way, so it seems preferable to build a model for the mass 
distribution that takes into account this important fact. 
Using this more appropriate model for Draco, we have 
obtained the 7-ray flux profiles for the case of a cusp and 
a core DM density profiles (both scenarios are equally 
valid according to the observations). To do that, we first 
estimated the best-fitting parameters for each density 
profile by adjusting the solutions of the Jeans equations 
to velocity moments obtained for the Draco stellar sam- 
ple cleaned by a rigorous method of interloper removal. 
Apart from this recommendable and useful update on 
the best DM model for Draco, one important conclusion 
can be extracted concerning the absolute 7-flux: for both 
cusp and core DM density profiles, the flux values that we 
obtain are very similar for the inner region of the dwarf, 



i.e. where we have the largest flux values and signal de- 
tection would be easier. 

There is, however, a way to distinguish between a core 
and a cusp DM density profile. The crucial points con- 
cerning this issue are the sensitivity and the PSF of the 
telescope. If the telescope resolution is good enough (and 
we reach the required sensitivity) a distinction between 
both cusp and core models may be possible thanks to 
the shape of the flux profile in each case. However, if the 
PSF of the instrument is poor, its effect could make it 
impossible to discriminate between different flux profiles, 
i.e. different models of the DM density profile. In any 
case, to be sure that the signal is due to neutralino an- 
nihilation, we will need to have a PSF good enough to 
be able to resolve the source (i.e. we will need to detect 
with a good resolution at least a portion of the flux pro- 
file large enough so we can conclude that it belongs to 
neutralino annihilation). This fact together with a char- 
acteristic spectrum are the unequivocal traces of a DM 
source. Both issues are of course totally valid not only for 
Draco but also for any other target, and therefore they 
should be always taken into account. 

We may think that for Draco the effect of the PSF 
is not especially important, since current experiments do 
not seem to reach even the required sensitivity to detect a 
7-ray signal coming from the center of the dSph. Indeed, 
the limiting factor here is the expected flux and not the 
PSF of the telescope. Nevertheless, we should account 
for the PSF of the instrument in any case, since we do 
not know a priori the expected gamma flux and we will 
not have a realistic prediction unless we include it in our 
calculations. Only then we will know the exact shape of 
the flux profile as could be observed by the telescope and 
therefore we can evaluate the real chances of detection. 
Furthermore, we must note that the inclusion of the PSF 
corresponds to a more general analysis not only valid 
for Draco but also for any target (the case the Milky 
Way for example where the role of the PSF is critical 
to understand the origin of the gamma emission in its 
center). 

Some estimations concerning flux profile detection 
prospects for the MAGIC and GLAST experiments have 
been also shown. According to these calculations, even in 
the most optimistic scenario (i.e. with the highest values 
of Jsusy allowed for the MSSM particle physics model 
adopted here) a 7-ray signal detection from Draco seems 
to be really hard at least for the current IACTs and for 
the GLAST satellite. In addition, it would be really dif- 
ficult to improve our expectations for Draco, since for 
example larger integration times will not improve drasti- 
cally the sensitivity lines of these instruments up to the 
level allowing successful detection (we need around five 
and three orders of magnitude more sensitivity for success 
with MAGIC and GLAST respectively). The uncertain- 
ties coming from astrophysics, as already mentioned, are 
negligible compared to those coming from the particle 
physics, so either by choosing a cusp or a core DM den- 
sity profile for Draco (or even modifying this for a steeper 
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one) we will not be able to increase the 7-ray flux up to 
the level where we can expect a signal detection. 

We also explored the prospects of an excess signal de- 
tection (i.e. we are not interested in the shape of the 
gamma ray flux profile, only in detectability) for MAGIC 
and GLAST, but we reached the same strongly negative 
conclusions as those obtained for the flux profile detec- 
tion. Both instruments are very far from the required 
sensitivity so it seems that they do not have any chance 
of successful detection. But even more, in the case of an 
(expected) unsuccessful detection we would not be able 
to put any useful constraints on the particle physics in- 
volved (the uncertainties from different dark matter den- 
sity profiles and other astrophysical considerations are 
negligible compared to this). None of the actual allowed 
values for fsusY could be rejected if we find no detection. 
This makes Draco even less attractive for those current 
gamma ray experiments that try to find DM traces. 

From these negative results it seems that, at least for 
current experiments, Draco does not represent a good 
target for DM searches (although some important en- 
hancements mechanisms in the 7-ray signal coming from 
neutralino annihilation in Draco have been proposed in 
[l|| ). If we want to find unequivocal traces of non- 
barionic dark matter in the Universe, an effort should 
be done in order to find other more promising DM tar- 
gets. In this context, it is interesting to note that the 
number of known dSph galaxies satellites of the Milky 
Way has increased considerably in the last few years, do- 
ing possible that some of them may offer more optimistic 
detection prospects than Draco itself [5l|. A quantita- 
tive study should be done in this direction. The search 
of DM subhalos in the solar neighborhood (see e.g. (52J) 
or exploring the IMBHs scenario [53[ may represent also 



other challenging possibilities, although these one more 
speculative for the moment. 

Finally, it is worth mentioning that IACTs that join a 
large field of view with a high sensitivity will represent 
the future in this field and will provide a next step in 
DM searches. In this context, the Cherenkov Telescope 
Array (CTA) project will be specially important, with 
a threshold energy well below 100 GeV, a wide spectral 
coverage up to 100 TeV and a factor of 5-10 more sensi- 
tivity than present experiments in the mentioned range. 
Also GAW, a R&D experiment under development, with 
an energy threshold ~ 700 GeV and a 24° x 24° field of 
view, will constitute another complementary attempt in 
the same direction, although under a different approach. 
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